The One-Dimensional KPZ Equation: an Exact Solution and its Universality 
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We report on the first exact solution of the KPZ equation in one dimension, with an initial 
condition which physically corresponds to the motion of a macroscopically curved height profile. 
The solution provides a determinantal formula for the probability distribution function of the height 
h(x,t) for all t > 0. In particular, we show that for large t, on the scale t 1 ^ 3 , the statistics is given 
by the Tracy- Widom distribution, known already from the theory of GUE random matrices. Our 
solution confirms that the KPZ equation describes the interface motion in the regime of weak driving 
force. Within this regime the KPZ equation details how the long time asymptotics is approached. 



The motion of interfaces persists as a fascinating topic 
of statistical mechanics. One particular, intensely stud- 
ied case are nonequilibrium growth processes, which are 
governed by local rules. In the seminal work [1] Kar- 
dar, Parisi, and Zhang proposed a model equation, now 
called KPZ equation, to investigate the dynamic scaling 
of such growing interfaces. KPZ argued that a grow- 
ing interface has a statistically self-similar structure with 
universal scaling exponents for the interface width and 
for the transverse correlation length. In particular, they 
predicted that a one-dimensional interface has fluctua- 
tions which grow as t 1 ^ 3 , in contrast to the i 1 / 4 statis- 
tical broadening of an equilibrium interface. Details of 
the KPZ scenario have been investigated through Monte 
Carlo simulations of simplified stochastic growth models, 
like Eden growth, random deposition, and polynuclear 
growth. Universal scaling exponents were confirmed, 
supporting that all such growth models constitute the 
KPZ universality class. On the theoretical side a variety 
of techniques have been devised, useful also in other areas 
of nonequilibrium statistical mechanics, as dynamic and 
functional renormalization group, mode-coupling theory, 
exact solutions, and more. For the development up to 
1995 we refer to the book [2] and the reviews [3l [4] . How- 
ever, so far the KPZ equation itself has resisted analytical 
handling. We will report here on the first exact solution. 
Physically it describes cluster growth in a thin film. The 
solution also provides a better understanding of the uni- 
versality of the KPZ equation. 

We will be concerned with the one-dimensional KPZ 
equation only. The interface location is then described by 
a height function h(x, t) depending on time t and location 
x on the real line. The KPZ equation reads 

The first term is the slope dependent growth velocity 
of strength A, A > for convenience. The Lapla- 
cian smoothens the interface with diffusivity v > 0. 
r\ is normalized space-time white noise with correlator 



(77(2, t)r)(x' , t')) — 8{x — x')5(t — if) and \jD is the noise 
strength. 77 models the random nucleations at the in- 
terface. We will choose initial conditions such that the 
height profile remains curved on the macroscopic scale. 

In the breakthrough contribution [5] Johansson firstly 
succeeded to compute a universal probability distribu- 
tion function (pdf), see also the related work [SI [7]. He 
studied a discrete growth model, known as single step, 
with wedge initial profile. Most surprisingly, he discov- 
ered that the pdf for the random amplitude of the height 
is the Tracy- Widom distribution [5] , first obtained in the 
context of the large N statistics of the largest eigenvalue 
of a GUE random matrix [§]. The Tracy- Widom pdf, 
Ptw(s), is defined through the determinant of a symmet- 
ric operator acting on functions on the real line, where, 
in principle, the determinant is defined by the product 
of the eigenvalues. More explicitly, we denote by P s the 
projection operator onto the interval [s, 00) and by K A i 
the Airy kernel 

/>oo 

K Ai (x,y)= / dwAi(x + w)Ai(y + w) (2) 
Jo 

with Ai the standard Airy function. Then 

Ptw(s) = -r- det(l - P s K Ai P s ) . (3) 
as 

One can show that tr[P s KAiP s ] < 00 |8j and hence the 
infinite product makes sense. 

In a spectacular recent experiment on electroconvec- 
tion [TU] the Tracy- Widom statistics is verified down to a 
scale of height samples with probability 10~ 4 . In this ex- 
periment a thin film of liquid crystal is electrically driven 
to a turbulent phase, the unstable DSM1 phase. One 
then plants through a laser pulse a sharply localized seed 
of topological-defect turbulence, the stable DSM2 phase. 
It grows isotropically over a time window of 30 sec to 
the maximal size of 1.6 cm in diameter. The experiment 
also investigates how the Tracy- Widom distribution is 
approached for long times. One observes that the cumu- 
lants 2,3, and 4 have already reached their Tracy- Widom 
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FIG. 1. (a) A typical realization of the droplet height func- 
tion, (b) A directed polymer configuration. 



value, while the average still slowly decays as i -1 / 3 to- 
wards its stationary value, an observation so far without 
theoretical explanation. 

Exact solution. — The curved height profile can locally 
be approximated by a parabola. In idealization, we there- 
fore choose our initial conditions such that the aver- 
age height is exactly parabolic, which can be achieved 
through the initial sharp wedge 



h(x,0) = -\x\/S, <J<1. 



(4) 



Since for short times the nonlinear term dominates, the 
solution h(x, t) spreads rapidly into the parabolic profile, 
h(x, t) ~ -x 2 /2Xt for \x\ < Xt/S and h(x, t) ~ (\/26 2 )t- 
\x\/S for |x| > Xt/S. Because of the noise term in ([I]) 
this profile has superimposed random fluctuations whose 
amplitude grows as i 1 / 3 . A typical realization is shown 
in Fig. 1(a). 

We consider the height statistics at one prescribed 
point x. Then, for every t > 0, 



(\/2v)h(x,t) = -x 2 /Avt 



±1? 

12 It 



2 log a + 7*6= (5) 



where 



7t = (a 4 i/t) 1/3 : 



(2 V )- 3 / 2 XD 1 / 2 . 



(6) 



The first three terms of ^ are deterministic, in partic- 
ular one notes the inverted parabola already mentioned. 
The logarithmic term reflects the scale invariance of the 
initial sharp wedge. All information on the fluctuations 
of h(x,t) is encoded in the random amplitude which 
extends the Tracy- Widom amplitude, £t\Vi to finite time 
t. The prefactor j t confirms that the fluctuations grow 
as t 1 / 3 . The pdf of £t is given by 



Pt(s) 



du lt e lt(s - u) exp [ -e 



7t(s-«)l 



x ( det(l - P u (B t - P Ai )P u ) - det(l - P u B t P u )) . (7) 

Here Paa has the integral kernel Ai(x)Ai(?/), which can 
be viewed as a one-dimensional unnormalized projection, 
and B t has the kernel 



B t (x,y) 



dw(l~c-~ /tW )- 1 Ai(x+w)Ai(y+w). (8) 



£ t is independent of x and depends on t only through 
the dimensionless parameter 7 t . In the limit t — > oo the 
Gumbel density in Q tends to 8(s — u) and B t tends to 
= Ka\- The difference of determinants in ^ is then 
precisely prw(s)- Hence £ t tends to £tw as t — > oo. In 
case of a macroscopically curved profile, it is expected 
that for any growth process in the KPZ class the pdf of 
the height fluctuations is Tracy- Widom in the long time 
limit. We have thus established that the KPZ equation 
is in the KPZ universality class as regards to the pdf of 
the height fluctuations. 

Our exact solution also provides an explanation for the 
observed i -1 / 3 relaxation to stationarity. The difference 
of determinants in is of order t~ 4/3 pi] . Hence the 
leading correction to /?tw( s ) is the shift due to the non- 
zero mean of the Gumbel distribution which is of order 

r 1 / 3 . 

To represent the solution to ([T]) in a more tractable 
form we introduce the Cole-Hopf transformation as 



Z(x,t) = exp[(\/2v)h(x,t)} 



(9) 



At the expense of multiplicative noise, Z(x, t) then sat- 
isfies the linear equation 



9 7 ^ 7 



{\VD/2v)r)Z . 



(10) 



with the normalized initial condition Z(x,0) 



lim <5 ^o((2^/A)25)- 1 exp[-(A/2i/)|x|/(5] = 5(x). (|10]) is 
solved through the Feynman path integral 

Z(x,t) = E (exp a f darj(b(s),s) 8{b(2vt) - x)j . 



(11) 

Here Eo(-) is the Wiener integral over all paths of an 
auxiliary Brownian motion b(t), starting at and with 
variance t. In principle, the solution of the KPZ equation 
with sharp wedge initial condition is defined by 



h{x,t) = (2v/\)\ogZ(x,t) . 



(12) 



However, as written (Z(x,t)) — oo because of ultraviolet 
divergencies. As in quantum field theory one thus has 
to introduce a suitable cutoff to be removed through a 
renormalization scheme. We briefly describe four distinct 
variants, all providing physical insight to the interpre- 
tation of the KPZ equation. But only the last variant 
carries the computational power to obtain the exact so- 
lution. 

(i) Colored noise. — The expression ( 11 ) is only formal, 



since the white noise rj is very rough and hence the action 
in ( |ll[ ) relies on an undefined integral. One way to im- 
prove the situation is to substitute in the KPZ equation 
@> the noise T](x,t) by a noise rj K (x,t), which is colored 
in x with width k -1 but still white in t, such that rj K — s- r/ 
as k — y oo. This modification induces the uniform trans- 
lation of the height profile by v K t with v K ~ k. Going 
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to the moving frame of reference, the limit k -> oo of 
the height profile is well defined [121 S3] and yields the 
proper interpretation of (11). 

(ii) Directed polymer in a random potential. — The in- 



been accomplished by Tracy and Widom [T7] [TB] . Also 
available is the corresponding result for the PNG droplet 



The Cole-Hopf solution (11) arises at the crossover 



tegral in (11) is discretized. Then the Brownian motion scale to weak asymmetry (WASEP). More specifically, 

one assumes q = \ + y/e(3, (3 > 0, e > 0, and e « 1. 



path b(t) is replaced by a directed polymer u), which for 
convenience is placed on the two-dimensional lattice Z 2 . 
It starts at 0, makes only up/right moves, and ends at 
the lattice site (tN,tN), see Fig. 1(b). Independently 
for each site (i, j) G Z 2 , there is a unit Gaussian random 
potential rj(i,j). The energy of the directed polymer is 
E(lj) = where the sum is over the potentials 

along the polymer lo. The discrete approximation to ( 11 ) 
reads then 

Z tN = J2 exp[-PE{u)]. (13) 

w:(0,0)~>(tJV,tJV) 

The partition function is point-to-point, since both end- 
points of the polymer are fixed. In the limit N — > oo, 
/3 -> with /3 4 A fixed, Z tN -> Z(0,t) and Z(0,t) as de- 
fined in (i) [TJ] . From this perspective the KPZ equation 
is the weak noise limit of the directed polymer, see [TS] 
for details. 

(Hi) Attractive 5-Bose gas |16j . — The n-th moment, 
(Z(x, t) n ), of the partition function pTj ) can be expressed 
through the propagator of n quantum particles on the line 
interacting through an attractive ^-potential. All n parti- 
cles start at and propagate to x at time t. Working out 
(Z n ) yields the potential —\a 2 Y^7,j=i $( Xi ~ x i)- Thus 
the renormalization corresponds to merely removing the 
self-energy through normal ordering. Thereby one arrives 
at the same random partition function as constructed in 
(i) and (ii). 

(iv) Single step growth model. — This is a stochastic 
evolution model for the integer valued height function 
h{j,t), j € Z, satisfying the single step constraint 



\h(j + i,t)-h(j,t)\ = i. 



(14) 



In our case the initial height profile is the wedge h(j, 0) = 
— A random sequential update rule is used: Indepen- 
dently the height at a local minimum is increased by 2 
with rate p and at a local maximum decreased by 2 with 
rate q, p + q = 1 to set the time scale, and q > p in 
our case corresponding to A > 0. The height differences 
are then governed by the partially asymmetric simple 
exclusion process (PASEP). This is a stochastic particle 
system on Z, where there is at most one particle per site. 
Particles jump with rate p to the right and rate q to the 
left under the constraint of the exclusion rule. The initial 
wedge corresponds to the 0-1 step initial particle config- 
uration, for which all sites to left of the origin are empty 
and to the right of the origin are filled. In [5] the Tracy- 
Widom fluctuations are proved in case p = 0, q = 1 (the 
TASEP), including a discrete time parallel random up- 
date rule. The intricate extension to0<p<q<l has 



The correspondingly adjusted time scale is order e~ 2 
and space scale is order e" 1 , while the height is of or- 
der e -1 / 2 . If h e (j,t) denotes the WASEP random height 
profile, then the single step partition function is defined 
through 



^stcp 



(j,t) =exp[(\og(p/qW(j,t)}. 



(15) 



As 
Z 



the correctly 
the limit e 



centered 
^ to 



shown in [T2J IT5] . 

stcp (e~ 1 x,e~ 2 t) converges in 
Z(x,t), constructed already before in (i) - (iii). 

We are now in a position to indicate, rather roughly, 
how our exact solution is obtained. A more detailed ex- 
position can be found in [H1[2T]. For the 0-1 step initial 
particle configuration, Tracy and Widom [17] recently 
provided a contour integration formula for the proba- 
bility distribution of the position of the m-th particle at 
time t valid for arbitrary q. In |21[ we use their formula as 
starting point and study the scaling limit of the PASEP 
as above, namely asymmetry y/ef3, space 0(e~ 1 ), and 
time (D(e~ 2 ), to arrive at a pdf given as an integral over 
the difference of two determinants as in (f7l), including 
a deterministic shift of order e~ x and a subleading logs 
correction. As discussed in the companion paper [llj . in 
addition one has to analyze the average (Z^ tcp (j 7 t)} for 
the WASEP, thereby to determine the appropriate cen- 
tering. One finds terms of order e^ 1 and loge, which 
are precisely cancelled by the corresponding terms ap- 
pearing in the WASEP scaling limit. Combining both 
results, one arrives at (|5]) - ((8)). The KPZ parameters 
are fixed through the WASEP as 2v = 1 and /3 = A. At 
x = the WASEP average density is 1/2 and hence the 
noise strength D = 1/4. By varying x one can tune the 
noise strength to p(l — p), where p is the local average 
density. In principle it suffices to consider x — and to 
deduce the parameter dependence from ( |11| . 

For a more detailed information on the pdf pt{s) one 
has to rely on a numerical evaluation of (J7|, for which 
purpose it is convenient to regard, e.g., P u B t P u as a large 
matrix. One then evaluates B t (x, y) at a judiciously cho- 
sen set of grid points Xj, j = 1, . . . ,n, u < Xj, with n 
of order 60 to 120 and thereby defines the n x n matrix 
Bij = B t (xi,Xj), i,j = 1, ...,n, for which det(l — B) 
is obtained by a standard routine. Since the kernel of 
^ is a smooth function, the computation is fast and the 
errors are small [25]. In Fig. 2(a) we display the pdf 
Pt{s) for j t = 2, 5, 10 obtained by the method described. 
For even smaller values of 7 t the computation becomes 
slow because of strong oscillations of the kernel B t (x, y) 
resulting from the singularity at w = of the defining 
integral 
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FIG. 2. (a) The pdf defined by @ for 7t = 2(-),5(+), 10(-) 
and the Tracy- Widom pdf defined by (|3| (full curve), (b) 
The TASEP height statistics at 10 3 MC steps and the Tracy- 
Widom pdf. 



Universal properties. — The exact solution has been ob- 
tained through a limit of vanishing asymmetry. More 
generally, one could consider some growth model in the 
KPZ class with a tunable asymmetry, denoted by ft. 
Physically the asymmetry is possibly small, but fixed. 
Time is measured in Monte Carlo time steps. For times 
of order j3~ 2 the nonlinearity plays no role, yet, and the 
height fluctuations are approximately Gaussian. On the 
time scale /3 -4 , one crosses over to the KPZ solution 
([5]) with non-Gaussian statistics, which asymptotically 
is well approximated by the Tracy- Widom distribution. 
Of course, a clean separation of time scales is expected 
only for (3 <C 1. 

To have a test case we performed MC simulations for 
the PASEP with strongest possible asymmetry, namely 
q = 1, p = Q (the TASEP). For the 0-1 step initial condi- 
tion we sampled the height at the origin, h q= i(0,t). The 
scale linear in t is given by —(q — p)t/2. For t = 10 3 
MC steps the distribution of h q= i(0,t) is concentrated 
approximately in an interval of 60 lattice units. As dis- 
played in Fig. 2(b), even on the discrete level the Tracy- 
Widom distribution is still an accurate approximation, 
which becomes almost indistinguishable upon shifting it 
by 0.13 units to the right. The t" 1 / 3 shift to the right 
was noted before numerically from the recursion relations 
of the PNG droplet [23] and is also observed in the ex- 
periment [10]. According to the exact solution, for large 
t, p t is shifted by 7 t _1 (— 0.577 + 21oga) relative to pxwi 
where —0.577 is the mean of the Gumbel density. Since 
a is proportional to A, this suggests that the prefactor 
of the t~ 1 / 3 shift for the PASEP should under go a sign 
change. Indeed, running the MC simulation for smaller 
values of q, one finds that pt is shifted to the left side 
of ptw f° r 1 0.78. The approach as i" 1 / 3 seems to be 
universal and is expected to hold for any growth model in 
the KPZ class. However the precise first order correction 
to ptw(s) will depend on the particular model. 

Conclusions. — We have computed the exact probabil- 
ity distribution function for the height h(x, t) of the KPZ 
equation with narrow wedge initial profile. The long time 
limit is the Tracy- Widom pdf, in accordance with the re- 
sults for discrete growth models. But in addition, based 
on our exact solution, we now better understand how the 



long time asymptotic is approached. For growth models 
with tunable asymmetry the KPZ equation is an accurate 
approximation for weak asymmetry. 

Note. — After submitting our letter, three distinct 
groups [T31 H21 [H] posted independent but related work 
on the arXiv. 
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knowledges the support from KAKENHI (9740044) and 
H.S. from Math-for-Industry of Kyushu University. 
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